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G^ '. Abstract 



Q^' We show explicitly that there is particle creation in a static spacetime. This 

0-) ■ 

is done by studying the field in a coordinate system based on a physical 
principle which has recently been proposed. There the field is quantized by 



^ 



X 

C^ I decomposing it into positive and negative frequency modes on a particular 

spacelike surface. This decomposition depends explicitly on the surface where 

the decomposition is performed, so that an observer who travels from one 

surface to another will observe particle production due to the different vacuum 

state. 
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I. INTRODUCTION 

Before one can quantize a free field propagating on a curved background one must have 
a unique means of splitting the field into positive and negative frequency modes. This 
procedure requires a unique definition of time with which to perform the decomposition. 

There have been many attempts to define the vacuum for a free field propagating in a 
nontrivial spacetime. It is a common feature of all these attempts that the choice of vacuum 
is determined by a particular choice of time coordinate. This is true even for such general 
quantization procedures as Deutsch and Najmi |T|] although there the dependence is not 
explicit. Instead they require a foliation of spacetime by a family of spacelike hypersurfaces 
which, in essence, defines "instants of time" and the normals to these surfaces define the 
"direction of time". The choice of time coordinate, in most computations, has usually 
been based on calculational convenience and not on a local physical principle. The physical 
principle which we use here is that suggested by Capri and Roy 0. The principle is that 
on the surface of instantaneity a 1+1 dimensional Poincare algebra including the Killing 
equations for the generators should be valid. 

In a globally hyperbolic spacetime with one timelike and one spacelike dimension the 
surface of instantaneity, in this case a line, for a given observer is given by the particular 
spacelike geodesic which passes through the point at which the observer is located and is 
normal to the observer's timelike worldline. The direction of time on this surface is then 
defined to be everywhere normal to this spacelike geodesic. It has been shown that this 
definition of time is the unique one to obey the physical principle given above 0. 

The coordinates on the spacelike surface are chosen, for convenience, to be Riemann coor- 
dinates based at the observer's position, although any other coordinatization of the spacelike 
surface will do. To define the direction of time Gaussian Geodesic Normal coordinates are 
then constructed on this spacelike surface so that the time coordinate of some point off the 
surface is just given by the distance from the point to the spacelike surface. 

When one expresses the metric in terms of these new coordinates one finds that the 



metric has the form 



ds'^ = dt^ + giidx^, (1.1) 

with gii < and 5^11 = — 1 + 0{t^) near the origin of coordinates. 

To now decompose the field into positive and negative frequency modes we impose ini- 
tial conditions that force the field to have the correct time dependence (exp(— icjt)) in the 
neighbourhood of this spacelike surface. To ensure this correct time dependence we impose 
the initial conditions, 

(/)+(t,x)|i=o =An(0,x) and (<9t(/)+(t, x)) |i=o = -«cJn(0)A„(0,a;). (1.2) 

where the v4„(0,x) are the spatial eigenmodes of the Laplace-Beltrami operator on the 
surface t = where the decomposition is to be performed. 

In section II we construct the metric in terms of these physically preferred coordinates 
and write out explicitly the boundary conditions which determine the positive frequency 
modes at the surface t = 0. In section III a complete set of modes for the Laplace-Beltrami 
operator for a massive scalar field is obtained and in section IV the orthogonality relation 
satisfied by these modes is calculated. In section V we then use these orthogonality relations 
to impose the physically relevant boundary conditions which were calculated in section II. 
In section VI the actual particle creation due to the presence of the gravitational field is 
calculated for a stationary observer with respect to this spacetime, in terms of the original 
static coordinates. 

II. THE PREFERRED COORDINATES 

The spacetime we are interested in is described by the metric 

ds" = a(X)dT^ - ^^ (2.1) 

a[Xj 

where 



a{X) = l-exp{-q{\X\-r)). (2.2) 

To simplify the technical discussion later on we choose r such that exp(gr) < 2, then 
— 1 < a{X) < 1. This spacetime was first investigated by Witten as a 1 + 1 dimensional 
eternal blackhole spacetime. 

To construct the preferred coordinates we must solve the geodesic equations for this 
spacetime. The first integrals of the geodesic equations yield 

where ei = ±1 and e = —1 for spacelike geodesies and e = 1 for timelike geodesies. The 
spacelike geodesic which is perpendicular to the timelike vector -^(1,0), (oq = Q^(-'^o)) and 
which can be treated as the tangent vector to the worldline of the observer at Po{Tq, Xq), is 
given by setting Cq = and e = — 1. We can therefore see on this surface Sq that To = Ti. 
Ti is the T coordinate of the point Pi{Ti,Xi) which is the point at which the geodesic from 
the general point P{T, X) intersects this spacelike surface orthogonally. The preferred time 
coordinate t is given by the distance along the timelike geodesic connecting Pi to the general 
point P{T,X) which is normal to the surface Sq at Pi. This timelike geodesic is given by 
( ^73|) with Cq = tti = a(Xi) and e = 1 so that 

t = r dX'^ = ['' dX' , '' (2.4) 

Jxi dX' Jx, ^a{Xi) - a(X') 

One can also calculate the change in the coordinate T along this geodesic, 



T - Ti = T - To = / dT= dX' , =■ (2.5) 

•^^1 ^^1 a{X')Ja{Xi)-a{X') 



These two equations allow us to express the metric in terms of the coordinates t and Xi. The 
preferred coordinate x on 5*0 is now constructed using a 2-bein of orthogonal basis vectors at 
Pq, eo(Po) and ei(Po). With eo(Po) given by ^=(1, 0) the tangent to the observer's worldline 
and p^ given by the tangent vector at Pq to the geodesic connecting Pq to Pi, the Riemann 
normal coordinates ?]" of Pi are given by 



sp'^ = r/"e^(Po) (2.6) 

where s is the distance along the geodesic Pq — -Pi- Using e^C/j^ = rjais (Minkowski metric), 
and the orthogonahty of p^ to eo(-Po) "we have 

V' = sp^e'^iPo) r^' = -sp^e'^{Po). (2.7) 

The surface Sq is just the surface rp = Q and the coordinate x is 

x = r,'' = -sp^eliPo) = r' dX' ^ (2.8) 

The preferred coordinates (t, x) are then given by solving the above integrals for X > 
or X < 0. After choosing ei = — 1 one obtains for X > 0, 



and for X < 0, 



-2 V-1 + e9(--+^i) tan-i(V-l + e9(-^+^i)) 62 
g 

2tanh-i(^^^^)62 
+ v^-^"""""' + To (2.9) 



2 
X = - {tanh"^(v/a7) - tanh"^(v/a^)} (2.10) 



_ 2e^^^^?^tan-i(V-l + eg(-^+^^)) 

V — I J 



T 



2 V-1 + e-(9('-+^i)) tan-i(V-l + e9(^-^i)) es 



? 



X 



2tanh-i(^^^S)62 

V-^+--""""-" +ro (2.12) 

Q 

2 
- {-tanh"^(v^) + tanh~^(v^)} (2.13) 



^^ -2Un-(^TO^ i (2.14) 

e 2 g 



where 



€2 = ±1. 



(2.15) 



From these coordinate transformations we can see that x and t both run from — cxd to +cxd 
and cover the region of the original space corresponding to \X\ > r, the region outside the 
horizon. The region inside the horizon is shrunk to a point. Furthermore, the region between 
the observer and the horizon (r < |X| < |Xo|) is covered twice. 
In terms of the coordinates (t, x) the metric is now 



where 



2 J„2 



ds = dt — (1 + tp{x) tan [tp(x)]) dx 



(2.16) 



p{x) = -sech [-B(a:)] 



and 



B{x) = tanh ^ [y^a^ 



xq 



(2.17) 



(2.18) 



We can see at this point that in this coordinate system, which does have a physical basis, 
the metric no longer appears static. 

In these new coordinates the Klein-Gordon equation for a massive scalar field is 



d^(P + I {dM\9\)) M + -^dx {^\g''d^ ct) + m^(t) = 0. 



(2.19) 

We now define instantaneous eigenfunctions An{t,x) of the spatial part of the Laplace- 
Beltrami operator, such that 

1 



dx ^J\g\g d^] +m 



Akit.x) =Uk{t)Ak{t,x) 



(2.20) 



The positive frequency solutions of ( |2.19|) are then defined as those which satisfy the 
initial conditions 



(f)t{t,x)\t=o = Ak{0,x) and {dt(l)t{t,x))\t=o = -iuJk{0)Ak{0,x). 



(2.21) 



These initial conditions ensure that the positive frequency part of the field has the desired 
time dependence at t = 0. These positive frequency solutions form a vector space which is 
made into a Hilbert space using the standard Klein-Gordon inner product. 
From the simple form of the metric at t = we see that 

k k 

Afc(0, x) = sm{2-B{x)) or cos{2-B{x)) and ul{0) = (A;^ + m^) (2.22) 

With the positive frequency solutions defined in this way we can then write out the 
quantized field as 

v[/i = j^" dk-^^ {<Ptit,x)ak + 4^^*it,x)al} , (2.23) 

where the subscript 1 of the field simply denotes the surface on which the positive frequency 
modes have been defined. In this expression we have written ujk{0) as Uk and we will 
continue this practice. Unfortunately (|2.16| ) is too complicated to obtain the general form 
of the modes in terms of the coordinates {t,x). This is, however, not really a problem as 
the point of this approach is to find out what boundary conditions should be imposed. It is 
therefore sufficient to solve the field equations in whatever coordinate system is convenient 
and then express these solutions in the preferred coordinate system to impose the boundary 
conditions. 

III. MODES OF FIELD EQUATION 

From the form of (|2.1|) we can see that in terms of the original coordinates (T, X) the 
field equations are separable. For this reason we solve for the modes in these coordinates 
and then express the solutions in terms of the preferred coordinates using the coordinate 
transformations given above ( |2.9| - |2.14]) . In terms of the coordinates (T, X) the Klein-Gordon 



operator has the form 
1 



, ^, dT<P - {dxa{X))dx(l) - a{X)dj,<f) + m^0 = (3.1) 



By assuming a T dependence for the field of the form exp{—iUpT) we obtain the following 
differential equation, 

dx{a{X)dx(p) + {^^ - m2)0 = (3.2) 

a[A ) 

To construct a self-adjoint extension for this operator we are required to construct solu- 
tions which vanish at the horizon where a{X) = 0. 

By making a change of variable to 2; = 1 — exp(— g(|X| — r)) = a{X) we obtain the 
following differential equation in terms of z, 



z{l - zy^"{z) + (1 - z){l - 2z)^'{z) + {—- /x^)^(^) = 0. (3.3) 



z 



where 

p' = 4 a^d /i2 = ^ (3.4) 

qZ gZ 

We are interested in constructing solutions outside the horizon so that |X| > r and z > 0. 
As mentioned above we also require that the solutions vanish at z = 0. 
The two independent solutions to ( p. 2] ) are 

^ip(^) = z"(l-^)'F(a,6,c,z) (3.5) 

where F(a, b, c, z) is an hypergeometric function and 

n = ip 



I = iyp"^ — fi"^ 

a = n + I 

b=n+l+l 

c=l + 2n (3.6) 

and 

^2p{z) = z''{l-zyF{a,b,c,z) (3.7) 



where 



n = —ip 
I = —iyp"^ — /i^ 
a = n + I 

c=l + 2n. (3.8) 



We can now finally write out the desired solution to ( p.3|) 

^{p,X) = [^2pi0)^ipiz) - ^,p{0)^2piz)] e(X). (3.9) 

The general solution to ( |3.1| ) can then be written, 

^(T, X)= dp {(A(p)^(p, X) exp{-iujpT) + B{p)^{p, X) exp{iujpT))} (3.10) 

J ^ 

We now impose the initial conditions ( p.22| ) which then give some physical meaning to 
the decomposition of this field into positive and negative frequency parts. To explicitly 
impose these initial conditions it is first useful to find the orthogonality relation satisfied by 
the modes \l/(p, X). 

IV. ORTHOGONALITY OF MODES 

To find the orthogonality relation satisfied by the mode \E'(p, X) we follow standard 
Sturm-Liouville approach and recall that the mode satisfies, 

z 



z{l - zfF"{p, 2) + (1 - z){l - 2z)F'{p, ^) + (- - /i')F(p, z) = 0. (4.1) 



We can also write out a similar equation which is satisfied by the modes F*{k, z). If one now 
multiplies the equation for F{p, z) by F*{k, z) and the equation for F*{k, z) by F{p, z) and 
looks at the difference of the two equations one can see that after integrating over the range 
2; = to 2; = 1 and integrating the two terms by parts once we are left with the relation. 



/' rf,^^:^Mj^lM = lini i2_i (F'(p, z)F*{k, z) - F{p, z)F'*{k, z)) (4.2) 

JO z[\ — z) 2-^1 (p"^ — A;^) 

:F'(p, z)F*(fc, z) - F(p, 2)F'*(fc, ^)) |,=o (4.3) 



z~^l (p2 _ /^2^ 



Because of the boundary conditions satisfied by F{p,z) and F*{k,z) the second term in 
this relation is identically zero. The first term, as we show, is proportional to a delta 
function. This shows that these modes are orthogonal. To see that this expression is indeed 
proportional to a delta function we first smear it with a smooth function of p and show 
that the result is proportional to that function evaluated at p = k. When one attempts to 
evaluate the limit in the first term one finds that all the various terms are proportional to 
a common factor which produces the delta function, this factor is 



lim ^ ^ ^^ ^ (4.4) 

z^i {p- k) ^ 

To proceed we introduce a regularization (1 — zY and write, 

F{p,z) = lim F{p,z){l-zy. (4.5) 

The integral we must evaluate is, 
fiP) ,. 



lim/ (ipZl^lim(l-2)2^i(l-z)-^(V^^^-vV^)_(l_^)*(v/^^^-vV^)l (4,g) 

e^oj_oo p — k z^^ y } 

This shows that there is no contribution to the integral from the regions where |jo — fc| > R. 
In these regions the pole at jo = A; is not realized so one may interchange the order in 
which the limits are performed. These contributions then go to zero as the limit 2 ^ 1 is 
performed. We are then left with the integral 

lim lim f'^^^dp^^ lim(l - zf' ((1 - ^)-^(v/^^^-V^^^) _ (1 _ zY^V^^^^V^^^) 

(4.7) 

It is now convenient to make a change of variable to the variable x where 

X = ^■'°(l-;)(P-'-) (4^8) 

10 



We next expand the integrand in powers of [p — k) and find that as i? ^ we are left with 
the smooth function f{p) evaluated at the pole multiplied by a function of k. 

Using the above analysis for the orthogonality relations in z one can then write the 
orthogonality relations in X, 



l\X\>r a[X) 

= 5{p-k)\F{k)\^ (4.9) 



where 



V2(fc - v/F"^7?)r(l - 2ik). 



nJkqsmh.{2'Ky/k'^ — ^^) 



X sinh(7r(A; - JP - ^2)) sinh(^(A; + JP - ^i"^)) (4.10) 



V2{k + yF^7?)r(l + 2ik). 



T^ykq sinh(27rv^A;2 — /i^) 



X sinh(7r(A; - ^P - /i^)) sinh(7r(A; + ^P - /i2)). (4.11) 

V. FREQUENCY DECOMPOSITION AND THE VACUUM 

We now decompose the field into positive and negative frequency parts by picking out 
the positive frequency part of the field as that which satisfies the initial conditions in the 
preferred coordinates. This allows us to extract the annihilation operator for this field and 
thus define the vacuum state for this field on this particular spacelike hypersurface (i.e. the 
appropriate hypersurface which passes through the point (To^-^o))- The physically relevant 
question is, of course, how this decomposition depends on the point {Tq,Xq) which could 
represent the position of an observer. If this decomposition depends on the position of the 
observer then at some different position presumably the observer would observe some sort 
of particle density due to the change in composition of the vacuum state. To see this we 

11 



must impose the initial conditions relevant to the quantization on this surface. Recall the 
initial conditions 

(j)t{t,x)\t=o = A{0,x) and {dt(f)t{t,x))\t=o = -iuk{0)Ak{0,x) (5.1) 

where 

k 1 

Ak(0, x) = sm(2-B(x)) and uM) = ik'^ + m^)^ (5.2) 

We can thus write the general form of the solution which satisfies these initial conditions for 
this particular mode k 

^fc(r, X)= dp {iAkip)^ip, X) expi-ioOpT) + Bk{p)^ip, X) expiiUpT))} (5.3) 

J fi 

where we now regard T,X and z as functions of {t,x). This can be easily done given the 

coordinate transformations of section II. Because the initial conditions are imposed at t = 

we need only be concerned with the form of this field and its derivative normal to t = 

for z{t = 0,x) in order to evaluate the expansion coefficients A{p,n),A*{p,n),B{p,n) and 

B*{p,n). 

By using the orthogonality relations calculated in the last section we can write out the 

initial condition equations, 

\F{k)f {Ak{k)exp{-iuJkTo) + Bk{k) exp (iukTo)) 

'^*(k,X)sm(2!^B(x)) 
dX \;/ (5.4) 



dX \,J,,QT^' (5-5) 



\F{k)f{Bk{k)exp{-iuJkTo) - Bkik) exp{iuJkTo)) 

_^*(A;,X)sin(2^S(x)) 

l\x\>r^^ a(X)(f)|i=o 
In taking the time derivative of ( |5.3| ) one does not pick up a ^ because at t = this is 
zero. Again in these expressions it can be seen that we are still regarding z as z{0, x) and 
X is the inverse of this function in the integral. We have now determined A{p,n),A*{p, 
B{p,n) and B*{p,n) and can therefore decompose the field explicitly in terms of posit 
and negative frequency modes on this surface, 



nj, 
ive 
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v[/i= rdk^{<Pt,it,x)a,ik)+<P^^>it,x)a,iky} (5.6) 



where the extra subscripts denote the surface on which the frequency decomposition has 
been performed and the modes 0^^ (t, x) and 0^i(t,x) are the ones constructed with the 
expansion coefficients which satisfy (|5.4|) and ( |5.5|) . One may now define the vacuum relevant 
to this field on the surface (t = 0) in the usual way, 

ai{k) |0i) = \/ k (5.7) 

where again the subscript denotes "when" this is the vacuum state for the field. To see 
whether particles are created by the gravitational field in this spacetime one must look 
closely at how this state depends on the surface chosen. So at this stage all one needs is 
the normal derivatives, with respect to the spacelike surface, of the field on the surface. In 
Section VI the full transformation equations will be required. 

VI. PARTICLE CREATION 

To see whether particles are created by the gravitational field in the spacetime one 
must look closely at how the field decomposition depends on the surface chosen (i.e. the 
position of the observer). To obtain the spectrum of particles created one must calculate 
the Bogolubov transformation between the different annihilation and creation operators and 
look at the mixing of positive and negative frequency parts. To calculate the Bogolubov 
transformation we can just match the field from two different quantizations on a common 
surface. The easiest way to do this is to propagate one field to the surface on which the 
second is quantized. We can therefore write 

^i{t,x) = ^2{0,x') and dt'^,{t,x) = dti^2{t',x'))y=o (6.1) 

where t is the proper distance between the two quantization surfaces. This distance will, in 
general, depend on where on the surface one calculates the distance. To make the calculation 



13 



simpler we take Xq = Xq so that the observer is stationary with respect to the original 
coordinates where the metric is static. 

Because of the simple form of the modes at t = one can calculate the Bogolubov 
coefficients and write an expression of the form 

«2(^) = J '^P ("(^' k)ai{p) + I3{p, k)ai{py'^ (6.2) 

The particle density experienced by an observer travelling from surface 1 to surface 2 is then 
given by 

\l3{p,k)f (6.3) 

In 1 + 1 dimensions l3{p, k) in general has the form, 



In this equation the factors ^ and |p are required because we are matching the field's 
normal derivative is done with respect to the second surface. 

To calculate an approximate form of (3, valid for short time intervals, we expand the 
integrand about t = 0. To O(t^) we obtain, 

/■oo sinf— y) 2/r 

Pip,k) = - / dy \L_ tanh^(y) smi-y)p\x')uj,iT;, - T,f (6.5) 

where p{x') is given by (|2.17|) and we have changed variables from x' toy = B{x'). In getting 



from (|6.4|) to (|0| ) the second term of ( |6.4|) doesn't contribute to the integral because it is odd 
in y. It should be restated that this is particle creation observed by an observer stationary 
with respect to the original static coordinates. This can now be rewritten as, 

47r y uj^iUp J-oo cosh [yj q 

Several comments are in order here. Firstly, (3{p, k) is clearly non-zero so that particles 
are produced in this short time interval 5t = T^ — Tq. 
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Secondly, our approximation clearly only holds for Uk < q since the expansion breaks 
down for UkSt > 1. This means that we can only crudely estimate the number of particles 



produced in the time St since an ultraviolet cutoff of A; = y/q"^ — m? is required. 

Putting all this together we see that the momentum density of particles labelled by k 
produced in the time interval t is: 



nt{k)= rdp\P{p,k)f 
Jo 

J2K+2 r—r ^nc, A^ ^2 



q^5v ./uJ^ f°^ dp -K^ 



4tt q Jo ^pp^^^Trr? (f 



(p + A;)(^ + 3) b-A:)((£^ + 3) 



sinh^^i2±J. ginh^!^£ — '- 

q q 



(6.7) 



Within the spirit of the approximation, the total number of particles created in the time 
St with Uk < q is: 



Nt = / nt{k)dk (6.8) 

Jo 

This integral is finite, of course. If the upper limit is allowed to go to oo then the integral 
diverges linearly. This does not mean that the Bogolubov transformation is not unitarily 
implementable. Our approximations simply break down and our results are inconclusive. 
The difficulty arises from the fact that there are two time scales namely Ti = - and T2 = —. 
For a fixed Uk it is possible to expand in j; where T is the smaller of Ti, T2. However, if Uk 
is unbounded no such expansion is possible. 

VII. CONCLUSIONS 

We have shown that although a spacetime may be static this may not preclude particle 
creation which is a time dependent phenomenon [Q as the gaussian coordinatization may 
not be static. The only metrics which always lead to static Gaussian coordinates are those 
which have been dubbed "ultrastatic" by Fulling ||^. We have shown explicitly in this simple 
1+1 dimensional case how the choice of which coordinates should be used leads to some 
interesting results. In particular, the coordinates which are chosen via a physical principle 
seem to suggest that although the spacetime may be manifestly static in one coordinate 
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system these may not be the coordinates that one should use to quantize a field propagating 
in the spacetime. 

Unfortunately the analysis to find out whether the Bogolubov transformation is unitarily 
implementable was inconclusive. This is due to the fact that the approximate form of /3 which 
was analysed was not valid for large k. 
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